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ABSTRACT 


A model of Micro Multi Bore Hollow Fiber catalytic membrane reactor for the esterification process is developed. The catalytic 
membrane having a functional group is used to prepare 2-ethyl hexanoate from the reaction of methanol and 2-ethyl hexanoic acid. 
The implanted catalyst into membrane pores is playing a significant role in removing one of the products from the reaction mixture, 
thus breaks the thermodynamic equilibrium barrier, and increases the reaction conversion. The multi bores membrane reactor 
behaves as a plug flow reactor, considering the concentration is varying along one of the space coordinates (length), it is described 
by ordinary differential equations with and without respect of the axial dispersion. The membrane configuration with the perfect 
microcapillary channels gives regular flow pattern with minimum resistance. Experimental results, of an earlier project, are used to 
verify and validate the model; this approach gives a reliable tool for design, scaling-up and optimization. The developed model was 


based 84 


ARTICLE 


solved numerically. It revealed that the considering of the axial dispersion model is a reliable assumption proved appropriate 


matching with the predicted and calculating reaction conversion. 


Keywords: Multibores, plug flow reactor, membrane reactor, modeling, axial dispersion 


1. INTRODUCTION 


Esterification is a reversible reaction, conducted in the presence of a catalyst to promote the forward direction, the acidic media 
either homogeneous (liquid) or heterogeneous (solid) acids are usually used (lachon et al., 2003; Westermann and Melin, 2009), 
moreover, high temperature and long residence time are necessary to activate the heterogeneous acid catalysts are preferable due 
to the ease of separation, preventing corrosion reaction to get an optimum conversion. Using heterogeneous catalysts as ion 
exchange and Nafion resin can overcome some of the homogeneous catalyst problems, and give high reproducibility (Buonomenna 
et al., 2010; Jeong et al., 2003). The conversion is achieved to exceed the equilibrium conversion by selective removal of one of the 
products (Le Chatelier principle) (Levenspiel, 1999). 

The Membrane Catalysis has great respect comparing with traditional catalysts regarding its easiness to eliminate one of the 
products (usually water), which enhance the break of the thermodynamic barrier, thus rapidly obtaining high reaction conversion 
(Thomas et al., 2010; Shelepova et al., 2011; Plazaola et al., 2019; Buonomenna et al., 2010), it operates continuously, hence, the 
continuous regime gives excellent product quality and quantity. 

Polymeric membranes with Micro Multi Bores Hollow Fiber (MMBHF) are a new configuration, which, could combine the 
characteristics of flat sheet and hollow fiber membranes; therefore, it can tolerate elevated pressures, and it has a small footprint 
with a large surface area as hollow fibers membranes (ElZanati et al., 2016; Yingbo et al., 2015; Peng, and Wang, 2012; Xiao-Hua et 
al., 2013). It is considered the appropriate candidate to act as catalyst carriers due to its exceptional mechanical properties, high 
surface area, high packing density, easy module fabrication, outstanding stability, and reduced fiber breakage (Castanheiro et al., 
2006; Klofutar et al., 2010). 

According to Le Chatelier's Principle, the reaction conversion can exceed the equilibrium conversion by adding one of the 
reactants or removing one of the products, in this context, the water molecules are absorbed by a functionalized the grafted sulfonic 
group (Levenspiel, 1999), inside the membrane pores, where each pore is considered as a micro-reactor. The great contact between 
reactants and catalytic sites in the membrane reaches to high catalytic activity and minimum resistance (Xiao-Hua Maa et al., 2013). 

The objective of this article is to develop a mathematical model of the (MMCF) membrane reactor, conductingan esterification 
reaction.This model is verified using experimental results generated from a project funded by Science and technological 
development fund, Ministry of Scientific Research, Cairo, Egypt (ElZanati, 2017). 


Reactor Design 

Multi Bore Hollow Fiber Membranes 

Thirteen Membrane fibers, each one has seven micro bores, were earlier developed (El Zanati, 2017), Figure (1) demonstrates the 
cross-section of the elaborated fibers. 


Membrane bundle 


Sampling port 


Cross-section of 
MBHF membrane 


Membrane Reactor Experimental Set-up 


Figure 1 The Seven Bores Hollow Fiber Membrane reactor 
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Reactor construction 


The used catalytic membrane reactor is composed of 13 grafted fibers of 7 bores and 20 cm length contained in a glass tube, has a 
sample opening every 5 cm to study the kinetics of the esterification process Figure (1). 


Grafting of the membrane fibers 

A designed spinneret of multi holes is used to cast the fibers, prepared from Poly Ether Sulfone. The functional sulfonic group was 
grafted into pores of the microbore hollow fiber membrane using styrene as the monomer. The C-S double bond created by the 
styrene monomer can provide linkages that undergo chain polymerization (Thomas et al., 2010; Kumar et al., 2006). Activation step 
of the membrane is characterized by implanting a sulfonic group of negative charge on each branch of styrene polymer inside the 
pores; in which, 0.5N H2SO,4is permeated throughout the prepared membrane at 298° K. 


2. MODEL DEVELOPMENT 


Physical description of the reaction system 

The MMBHF catalytic membrane reactor is considered a multi micro plug flow reactors. The developed model relies on the 

geometry of the membrane, where, the lateral convective diffusion through the membrane wall occurs on internal pores area, which 

is measured by Brunauer-Emmet-Teller area, meanwhile, the axial dispersion may be affected the reaction propagation Figure (2). 
Fick's Law describes the axial dispersion flow through the perfect membrane cylindrical bores with the same radius. 

N; = —D, zal (1) 


Reactants FA + Methanol 


Lateral 
diffusion 


Axial 
Dispersion 


itn; 


Figure 2 The geometry of MMBHF catalytic membrane reactor and flow pattern 


Model assumption 

The following assumptions are adopted to develop the model. 

1- The MMBHF membrane reactor can be considered a plug flow reactor of radial homogeneities till exit, the fluid moves through 
the channel with uniform velocity. 

2- The flow inside the micro bores is considered perfectly mixed. 

3- In the plug flow, the convective flow is dominating besides of the dispersion of mass and heat-driven concentration and 
temperature gradients. 

4- Plug flow reactor is a steady-state reactor. 

5- Constant density. 
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6- Single reaction occurring in the bulk of the fluid (homogeneously), one of the products (water) is linked and adsorbed on the 


active functional sulfonic group, which was grafted inside the pores of the membrane walls. 

7- The reaction is isothermal and isobaric. 

8- The perfectness of the cylindrical micro bores is satisfactory to neglect the resistance. 

9- Mass transport of the reactants is axial and lateral-dispersed flow. 

10- Diffusion coefficients independent on the concentration, that is constant throughout the membrane thickness. 

11- The magnitude of dispersion is sufficient to account for the joint effect of all dispersive phenomena. The model is known as the 
dispersion model (T. Westermann et al., 2009). 

12- The active sites, in the membrane pores (SO3- groups), have an affinity only towards the water, the removal of water molecules 
shifts the esterification reaction to the forward direction according to the Le Chatelier principle. 

13- The unreacted alcohol, 2-ethyl hexanoic acid, and ester are not adsorbed on the membrane. 

14-The driving force of water transfer from the pore bulk to the pore surface was the difference between their concentrations at 
equilibrium. 

15- Fick's law describes the axial dispersion of mass flow (equation 1). 


Model derivation considering the axial dispersion 
Figure (2) demonstrates a longitudinal capillary section of the membrane and the grafted functional sulfonic group. The equation (2) 


describes the reaction. 

ky 
a+b 7? e+w (2) 

k, 
- The reaction rate constant can be written as (equation 3): 
r= k,C,Cy — k2C,.Cy (mol/cm2.s) (3) 
The forward reaction rate constant k; is a function of the reaction temperature; it is determined using the Arrhenius equation: 

-E 

k, =k,e /RT (4) 
The backward reaction rate constant kz is equal to ki/K. The equilibrium constant (K) is calculated using the following equations 


(Ichikawa et al., 2010): 


-AG 
InKzog = a (5) 
k= Kye “t/t (6) 


The values of AG and AH Table (1) are the difference of Gibbs free energy of products and reactants, and the difference of the heat 


of formation of products and reactants respectively. 


Table 1Gibbs free energy and the heat of formation of products and reactants 


Item Methanol 2-ethylhexanoic acid Water Meth. Eth. Hexanoate 
AG, kJ/mol -166.4 =25)\|1/ -228.61 =i | 
AH, kJ/mol = -238.4 -333.0 -285.8 -399.3 


Component mass balance: 


nip-+ OFS+ ykyS (Ci — kgCis) = mio ++ Ae (7) 
Considering that: 

nj = qCi and Nj = — io 

Cig + OFS + yiky S (Cj — kgCis) = qoCi + — — A,Dj oe (8) 


By rearrangement and dimensional correction: 
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dc; dA, dC 

giCie + 11S + y,ky (Ci — KaCis) = doCi + (S) vp St x ($5) GE — dh ADI GE (9) 
S = Sy — Ag, The internal pore surface area (cm?) at any time t. (10) 
ds = —dA, (11) 
Vp = Voo— Asx (12) 
dV, = —x dA, (13) 
Water adsorption, Nys: 
The change of water adsorbed on pore surface with channel length is: 
dnws d(CwsAs Ac 
“ss = _ = =e Vy Kiw (=) (Cy ~ KqgCws)S (14) 

Ac 
Gee dl < | a = Kw GC) Cw ~ KgCws)S (1 5) 


The change of the surface area covered with water (cm?) with length of channelis equal the change of water volume inside the pores 
divided by the water layer thickness; 


Ac 
dA, i Kw Gp Cw kaCws)S 


dl X Py (1 6) 
The change of water adsorbed on membrane surface with capillary length: 
Kw 2)Cws(Cw—KkaCws)S 

dCws _ 1 Ac = = i af 
Seen kw (32) (Cw — kaCws)S — | (17) 

Kiw(Cw-k Cws)S 
do = gr — MvGartasns (18) 
Rewriting equation (9) for each component (a, b, e and w) gives: 
dCa _ 1 = cha dA, 
a= (iy, [act — 1S — qoCa +x(H) Ca FE + dyADa Se am (19) 
dCy _ 1 = dA, d?Cy 
a = yy, [arcs rS — doCp + x(# 5) Cp al + dpAcDp ae | (20) 
dCe_ dA, 
<t = a [ace + 1S — qoCe +x(S) Co FE + dyADe Se <ce (21) 
dCw _ lw: 
a= Gy Tay [aeCwe + FS — Ki S (Cy — KaCws) — doCw +% (Se) Cw GE wt] (22) 


Ignoring the effect of axial dispersion 
A reforming of the previously developed model considering the ignoring of the axial dispersion effect on the change of reaction 
concentration along the length of the reactor. Equation (20) demonstrates the simplified model: 


gpCig + Oj(Sp — Ag) (Ki CC, — k2CeCw) + y;ki(So — As) (Cj — ka Cis) = Go (23) 


io i Ta 


3. RESULTS AND DISCUSSION 


The nonlinear model of esterification reaction accomplished in micro reactor of multi bore hollow fiber membrane is presented. The 
reactants are entered inside the micro bore hollow fiber lumen, and then dispersed axially. 
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Model solution 


The axial dispersion plug flow model for micro multi bores membrane reactor is developed assuming that each bore is functioning 
as a plug flow reactor. The model was solved for the change reaction conversion versus the independent spatial longitudinal variable 
(length, 1) at specific simulated state variables Cig (Cag,Cpo,Ceo, ANd Cwo), Cws, do, SoaNdAg, using Simulink, Matlab 52014. The earlier 
experimental results produced and published from the project titled “Novel Membranes with Multiple Hollow Holes (MMHH) for 
Reaction Catalysis to produce biofuel”, ID 9185, funded by Science and Technology Development Fund, Ministry of Scientific 
Research, Cairo, Egypt (El Zanati, 2017). 

The activation energy (E) for the esterification reaction of methanol with 2-ethyl hexanoic acid has been reported in the range 
from 42 to 53 7 kJ mol-'(Moya-Leon et al., 2006); accordingly, a value of 50 kJ mol" (average) is appropriate and used to estimate 
the reaction constants. 


Design parameters 

- Membrane bore diameter: 0.02 cm (ElZanati et al., 2015). 

- Total internal pore surface area, S. (of 13 fibers of 7 bores), 114.296 cm?. 

- Membrane pore volume, Vo; 0.57148 cm?. 

- The cross-sectional surface area; 0.028574 cm2. 

- qr 0.2105 cm? s~"; reaction temperature, 25°C (298 K); catalystloading capacity based on the IEC corresponding to astyrene mass 
fraction 0.12 during grafting formation [18]. 

- MR(nat/Npe), 1:1. 


Model verification 

The governing equations (9-22) were solved and verified using the experimental results obtained from the mentioned previous work 
(El Zanati et al., 2015). The sum of residual squares errors (SRS) (or Zerror 2 for N readings) was used to achieve the best fitting of 
results, Figure (3). 


Experimental result for A=xxx6e-10 


—G- Experimental reading 


——— Model Output, considering the axial dispersion 


oe 
oO 


Reaction conversion, X 
o o 
w n 


10 
Length,cm 

Figure 3 Curve fitting of reaction conversion calculated experimentally and that predicted using the model concerning the axial 

dispersion 


The results of the axial dispersion model were reasonably compatible with the experimental data simulated at the same reaction 
conditions at the model parameters gives the best are summarized as follows: 
- The initial esterification forward reaction constant, k (frequency factor), is equal to 2.92E-5 m?. mol™?.s7?. 
- The mass transfer coefficient of liquid water, ky: 1.764E-7 m.s7?. 
- The adsorption coefficient of water, kg: 6E-8 m7?. 
- The water adsorption thickness, X, was 5.8E-12 m. 

The maximum predicted reaction conversion X,obtained is 0.97606 with minimum error of 0.02 at equal 2.92E-5 as 


demonstrated in Figure (4). 
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0.0222 . - . ; : A . : 


|—Error Between Model and Experimental 
0.0222 |- 


0.0222 |- 


0.0222 - 


0.0222 |. 


Error Value 


0.0222 |- 


0.0222 |- 


0.0222 - 


A 


2.915 2.92 2.925 2.93 2.935 2.94 2.945 2.95 2.955 
ko «10% 


4 4 


0.0222 + 4 


Figure 4 minimum errors for best fitting of all parameters 


The solution of the second model approach (ignoring the axial dispersion) simulated at the same model parameters; does not fit 
the experimental results; where the maximum X,is 0.796 (18.45%) decrease compared with that obtained when axial dispersion is 
considered. Therefore, the supposed axial dispersion in the plug flow reactor is excellent approach, Figure (5). It is, also, noticed from 
Figure (3) that the reaction conversion is approached more than 90% at reactor height 5cm, then X, gives an asymptotic value; 
therefore, the length of the reactor can be designed in this case; shorter than 20cm, this behavior might be attributed to the large 
inner surface area of the membrane which reveals high number of pores (contain the active sulfonic group) in the micro-multi 
channels hollow fiber membrane. 


& Xexp ——xX,withdispersson ——xX, without disperssion 


o 
co 


o 
roa) 


o 
Bb 


o 
N 


Reaction conversion, X 


oO 


0 0.05 0.1 0.15 0.2 
Reactor length, m 


Figure 5 Conversion versus channel length for the two models 


The developed model with axial dispersion is clearly showed that it was valid and could be applied at other design/operating 
conditions with reasonable reliability and accuracy. The model is therefore suitable for design, replication, scaling-up, optimization, 
and parametric investigation. 

The model considering the axial dispersion is used to check and verify the concentrations of methanol and 2-ethylhexanoate 
(ester) concentrations with length (C, and C, ), Figure (6). 
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@ Ce,exp @ Ca,exp ——Ce,pred ——Ca, pred 
5000 
4000 
3000 
2000 


1000 


Concentration, mol/m3 


0 


0 0.05 0.1 0.15 0.2 
Length of reactor, m 


Figure 6 Best fitting of concentration of experimental results and the predicted by the model concerning the axial dispersion (Ca/Cp 
and Ce versus length) 


The verified model was also used to detect the change of the water concentration in product stream with length and compared 
with that obtained from experimental data. It is clear that the most of produced water is adsorbed on the membrane pores, while a 
trace percentage (7.27E-9%) of water released in output stream,Figure (7). 


A Cof water formed (adsorbed plus permeated 
-——= Cw permeated,pred) 


Concentration, mol/m3 


0 0.05 0.1 0.15 0.2 
Length of the reactor, m 


Figure 7 Change of concentration of Water in Permeate and concentration of water adsorbed on membrane 


4. CONCLUSION 


A plug flow with axial dispersion model was developed for esterification reaction of methanol and 2-ethylhexanoic acid in micro- 
multi channels hollow fiber catalytic membrane reactor. Another approach of ignoring the axial dispersion was developed. The two 
developed models are solved and verified using published experimental results. The initial reaction constant was adjusted to obtain 
good agreement between the reaction conversions determined from the model and experimental results. The simulation process is 
executed at the reaction conditions: MR of 1:1 (Methanol: 2-ethylhexanoic acid), T of 298°K, and IEC corresponding to 0.12 styrene 
mass fraction, and So 114.296 cm2. The model solution revealed that the consideration of axial dispersion in the micro channels plug 
flow reactor is convenient and satisfactory. The model will be used to study scale-up of the reactor to semi-commercial design size, 
and accomplish a preliminary feasibility. 
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Nomenclatures 


a : Alcohol (methyl) 

b : Acid (2-ethylhexanoic acid.) 

e : Ester of 2-ethylhexanoic acid 

w : Water 

i : Reaction components (alcohol, acid, ester and water), 

Ci : Bulk concentration of component i,molcm? 

Cis : Surface concentration of component i on pores surface, molm? 
Cus : Concentration of adsorbed water molecules on pore surface, molcm~ 
D; : Diffusivity of component i, cms"! 

E : Activation energy, kJmol"' 

AG : Difference of Gibbs free energy of products and reactants, kJmol"! 
AH : Difference of heat of formation of products and reactants, kJ mol" 
ky : Forward reaction constant, cm4 mol-'s~! 

ky : Backward reaction constant, cm4 mol-'s~1 

kl : Mass transfer coefficient of water in the bulk solution, cms”! 

kd : Adsorption coefficient of water inside the pores, cm" 

K : Equilibrium constant, dimensionless 

| : Length of channels, cm 

Ni : Mass flux of component, molcm* s7! 

nj : Number of components i, mol 

Nws : Number of adsorbed water molecules on pore surface, mol 


: Rate of reaction, mols“'cm-2 
on :-1 for reactants and +1 for products. 
:0 for i: a, b and e, while it equal —1 for i: w. 


qi : Flow rate, cm3s"! 

Pn : Molar water density, gmolcm™. 

dr : Feed flow rate, cm3s"! 

do : Output flow rate, cms"! 

SRS : Sum of residual squares 

R : Reaction gas constant: 0.008314, kJmol"! K"! 

T : Temperature degree, K 

Xp : Predicted reaction conversion 
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